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Abstract 

Riemann zeta function is an important object of number theory. It was also used for descrip- 
tion of disordered systems in statistical mechanics. We show that Riemann zeta function 
is also useful for the description of integrable model. We study XXX Heisenberg spin 1/2 
anti-ferromagnet. We evaluate a probability of formation of a ferromagnetic string in the 
anti- ferromagnetic ground state in thermodynamics limit. We prove that for short strings 
the probability can be expressed in terms of Riemann zeta function with odd arguments. 
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1. Introduction 



Riemann zeta function for Re(s) > 1 can be defined as follows: 

OO 1 

c(*) = E^ (i-i) 

n=l U 

It also can be represented as a product with respect to all prime numbers p 

as) = u^-p' s r 1 (i-2) 

p 

It can be analytically continued in the whole complex plane of s. It has only one pole, at s = 1 
and it has 'trivial' zeros at s = —2n ( n > 1 is an integer). The famous Riemann hypothesis 
[19] states that nontrivial zeros belong to the straight line Re(s) = 1/2. Riemann zeta function 
is useful for study of distribution of prime numbers on the real axis [18]. The values of Riemann 
zeta function at special points were studied in [21] , [22] . At even values of its argument zeta 
function can be expressed in terms of powers of 7r. The values of Riemann zeta function at odd 
arguments provide infinitely many different irrational numbers [20] . Riemann zeta function plays 
an important role, not only in pure mathematics but also theoretical physics. Some Feynman 
diagrams in quantum field theory can be expressed in terms of C( n )j see > f° r example, [1] . It 
appears also in string theory [2] . In statistical mechanics Riemann zeta function was used for the 
description of chaotic systems. This is large field with many publications. Important contributions 
to this field were made by Berry, Connes, Julia, Kac, Keating, Knauf, Odlyzko, Pitkanen, Polya, 
Ruelle, Sarnak and Zagier. One can find more information and citation on the following web cite 
http://www.maths.ex.ac.uk/ mwatkins/ . 

We argue that ((n) is also important for exactly solvable models. One of the most famous 
integrable models is the Heisenberg XXX spin chain. This model was first suggested by Heisenberg 
[3] in 1928 and solved by Bethe [4] in 1931 . Since that time it found multiple applications in solid 
state physics and statistical mechanics. Recently the XXX spin chain was used for study of the 
entanglement in quantum computations [6]. 

The Hamiltonian of the XXX spin chain can be written like this 

N 

h=y: + ofriVi + «i - 1 ) (i-3) 

i=i 

Here N is the length of the lattice and of, o\ , of are Pauli matrices. We consider thermodynamics 
limit , when N goes to infinity. The sign in front of the Hamiltonian indicates that we are 
considering the anti-ferromagnetic case. We consider periodic boundary conditions. Notice that 
this Hamiltonian annihilates the ferromagnetic state [ all spins up]. 

The construction of the anti-ferromagnetic ground state wave function \AFM > can be credited 
to Hulthen [5]. An important correlation function was defined in [12]. It was called the emptiness 
formation probability 

n 

P(n) =< AFM\ Yl Pj\AFM > 
i=i 

where Pj = (1 + <rJ)/2 is a projector on the state with spin up in jih lattice site. Averaging is over 
the anti-ferromagnetic ground state. It describes the probability of formation of a ferromagnetic 
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string of the length n in the anti- ferromagnetic background \AFM > . In this paper we shall first 
study short strings ( n is small), in the end we shall discuss long distance asymptotics ( at finite 
temperature). The four first values of the emptiness-formation probability look as follows: 

P(l) = \ = 0.5, (1.4) 

P(2) = -(1 - ln2) = 0.102284273, (1.5) 
3 

P(3) = \ - ln2 + -C(3) = 0.007624158, (1.6) 
4 8 

P(4) = - - 21n2 + — C(3) - — C(3) In 2 - — C 2 (3) 
55 85 

- 2|C(5) + 2|C(5) In 2 = 0.000206270 (1.7) 

Let us comment. The value of P(l) is evident from the symmetry, P(2) can be extracted from 
the explicit expression of the ground state energy [5]. P(3) can be extracted from the results of 
M.Takahashi [13] on the calculation of the nearest neighbor correlation. It was confirmed in paper 
[14]. One should also mention independent calculation of -P(3) in [15]. One can express -P(3) in 
terms of next to the nearest neighbor correlation 

< SfSfa >= 2P(3) - 2P(2) + l -P{l) (1.8) 

The calculation of P(3) and P(4) is discussed in this paper. 
The expression above for P(4) is our main result here. 

The plan of the paper is as follows. In the next section we discuss some main steps of the 
calculation of P(3) and P(4). The thermodynamics of P(n) for the non-zero temperature is briefly 
discussed in section 3. Then we summarize the results in the conclusion. 



2. General discussion of the calculation of P(3) and P(4) 

There are several different approaches to investigate P(n): 

• representation of correlation functions as determinants of Predholm integral op- 
erators described in detail in the book [10] 

• the vertex operator approach developed by the RIMS group [11] 

One can also mention the application of connection with other correlation functions, for in- 
stance, the correlation function < AFM\SfSf +n \AFM >. 

We shall use the integral representation obtained by Korepin, Izergin, Essler and Uglov [12] in 
framework of the vertex operator approach at the zero magnetic field: 

p( n )= f dXl f dX2 ,,,/ dXn TT(l| 1 y-g( y TT sinh7r(A fc - \j) 

Jc 2iri\i Jc 2m\ 2 " ' Jc 2iri\ n A a sinh ir\ a x <J-} k ^ n vr(A fe - Xj - i) ' 

The contour C in each integral goes parallel to the real axis with the imaginary part between 
and —i. 
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Recently such formula was generalized by de Gier and Korepin in paper [16] to the case, where 
averaging is done over arbitrary Bethe state [ with no strings ] instead of anti-ferromagnetic state. 

Let us describe in general the strategy we used in order to come to the answers (1.6) and (1.7). 
The integral formula (2.1) can be easily represented as follows: 

n r d\ 

p (n) = UJ c ^r U(X 1 , . . . , A n ) T(A l5 . . . , A n ) (2.2) 

where 

TT (\ \\ rii<fc<j<nSinh7r(A j - A fc ) 

^'■■■'^ - n ~ n^ =1 sinh"vrA, (2 ' 3) 

and 

T(Xl ,..., K) = (, 4) 

lh<k<j<n\ A j ~ A k — 1) 

As appeared we can make a lot of simplifications without taking integrals but using some simple 
observations. First of all, let us note that the function C/(Ai, . . . , A n ) is antisymmetric in respect 
to transposition of any pair of integration variables, say, Aj and A*,. This simple observation turns 
out to be very useful because 



n r 

j=i JC 



<1Xj U(X 1 ,...,X n )S(X 1 ,...,X n ) = (2.5) 



\Jc 2m 

if the function S is symmetric for at least one pair of A-s. 

The next observation is also trivial, namely, we can try to reduce the power of denominator in 
(2.4) using simple algebraic relations like 

(2.6) 



x(x + a) ax a(x + a) 

Combining these two simple observations one can reduce integration functions for P(3) to a 
sum of terms with denominators of power 2 and for P(4) to a more complicated sum of terms 
with denominators of power not higher than 3. 

In order to calculate the integrals one can close the contours in the complex plane by the infinite 
semi-circles either in upper half-plane or in the lower half-plane not changing the integrals. Then 
it is possible to apply Cauchy theorem using the following formulae 



dz f(z) 1 Id 
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= -^^--2^2) _/(*' + *) ( 2 - 7 ) 

f(il + e) (2.8) 



Ci 2ni sinh 3 nz 2n n 2 de 

dz f(z) 2 d Id 3 



Ci 2ni sinh nz 3n 2 de Air 2 de e 

for the cases n = 3 and n = 4 respectively where C\ is a small contour surrounding the point i I 
with an integer I in anti-clockwise direction. 

Then the integrals can be expressed in terms of the differential operator acting on some func- 
tions. For instance, for the case n = 3 
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where D is the differential operator 



7T 3 Id 2 Id 2 1 d 2 

8 7T dd\ 7T 2 de 2 7T 2 de 3 €ue2ie3 ^o 

sinh7r(e 2 — ei) sinh7r(e 3 — ei) sinh7r(e 3 — £2) (2-10) 

and 

00 00 00 

F( ei ,e 2 ,e 3 ) = El- 1 )' 1 El" 1 )' 2 E(-l) ,3 F(iZi + e u il 2 + e 2 ,tl 3 + e 3 ) (2.11) 

l 1= l 2 =0 l 3 =0 

Here all three contours were closed in the upper half-plane but in real calculations it turns out to be 
more convenient to close some of them in another direction taking into consideration appearance 
of an additional sign. 

It is not difficult to get generalization of these formulae to the case n — 4. So the problem is 
reduced to the calculation of sums like (2.11), expanding the result into the series in powers of 
e-s and applying the differential operator D. This procedure is straightforward but can be rather 
tedious especially for the case n — 4. Proceeding in this way we can come to the results (1.6) and 

(1-7). 

Let us note that both of these final answers appeared to be expressed in terms of the logarithmic 
function and the Riemann zeta function of odd arguments and do not depend on polylogarithms in 
spite of the fact that polylogarithm Li 4 (l/2) appeared in the intermediate stage of calculation. All 
coefficients before those functions in (1.4-1.7) are rational. Also they do not contain any powers 
of 7r which could be considered as Riemann zeta functions of even arguments. 

Our conjecture is that the final answer for any P{n) will also be expressed in 
terms of logarithm In 2 and Riemann zeta functions ((k) with odd integers k and with 
rational coefficients. 



3. Thermodynamics of P(n) 

If we had the exact answer for P{n) for any n we could calculate an asymptotics of P{n) when n 
tends to infinity. Unfortunately, for a moment we can not do this because we have P(n) only for 
n = 1,2,3,4. Nevertheless we can discuss a possible behavior of P(n) with n — > 00 using some 
other arguments. 

For non-zero temperature one can conclude that the asymptotics of the partition function in 
thermodynamic limit is as follows 

Z = < e^r > ~ e fr (3.1) 

where / is the free energy per site and is the length of the chain, it was evaluated in [7], [8] 
and [9]. In fact, for P(n) the n neighboring spins are frozen. Therefore one has the asymptotics 
of P(n) when n tends to infinity 

(l+<7?) H_ (N-u)f 
, \ < 117-1 efcT > e kT n f 

P{n) = I ^— = e-v (3.2) 

For zero temperature we expect Gaussian decay. 
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4. Conclusion 



We think that our work provide a link between integrable models and chaotic models. The same 
mathematical apparatus appears in the description of both kind of models. 

Let us repeat that the main result of this paper is the calculation of P(3) and P(4) (1.6- 
1.7) by means of the multi-integral representation (2.1). The fact that only the logarithm In 2 
and Riemann zeta function with odd arguments participate in the answers for P(l), . . . , P(4) 
and with rational coefficients before these functions allows us to suppose that this is the general 
property of P(n). One could compare the calculation of P(n) with the many- loop calculation of 
the self-energy diagrams in the renormalizable quantum field theory which can also be expressed 
in terms of ( functions of odd arguments [1] . 

Unfortunately, so far we have not got even a conjecture for P(n) but we believe that it is not 
an unsolvable problem. May be already after calculation of P(5) one could guess the right formula 
for a generic case P(n). It would give an answer to the question discussed in the previous section, 
namely, the question about the law of decay of P(n) when n tends to infinity. 

Also it would be interesting to generalize above results to the XXZ spin chain. Some interesting 
conjectures were recently invented by Razumov and Stroganov [17] for the special case of the XXZ 
model with A = —1/2. These conjectures would be supported if it were possible to get P(n) from 
the general integral representation obtained by the RIMS group [11]. 
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